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;^ ". The vectorial velocity is given as a function of the position of a particle 
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O . in orbit when a Newtonian central force is supplemented by an inverse 

p ,; cubic force as in Newton's theorem on revolving orbits. 
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approximate constant for orbits in general central potentials. A related 



Such expressions are useful in fitting orbits to radial velocities of or- 
bital streams. The Hamilton-Laplace-Runge-Lenz eccentricity vector is 

> . . 

oo ' generalised to give a constant of the motion for these systems and an 

O. 

^ ; vector is found for Hooke's centred elhpse. 

Qh! Introduction 

O , . . 

^ ! In a central orbit a particle with position (r, (p) has r'^dcjy/dt = h constant. Newton-*^ 

^ ■ 

^ ' pointed out that if = (f)/n, with n any constant, then r'^dcp/dt = h is also constant. 



He then enquired what extra radial force would be needed to make a new orbit with 
replacing but with the r{t) unchanged. Evidently h = h/n so the old radial equation 
of motion d'^r/dt^ — fi^r^^ = F will change to d^r/dt^ — h^r~^ = F. 

Thus the extra radial force required is F — F = — (ra~^ — l)h'^/r^ which is outward 
or inward according as n is greater or less than one. If the new motion were observed 
from axes that rotate (non- uniformly) at the rate Q = {n^^ — 1)0 then one would see the 
particle perform the original orbit unchanged. Thus from fixed axes the whole orbit may 
be thought of as revolving with angular velocity {w^ — 1)0 = (1 — n)0. This is Newton's 
theorem on revolving orbits. Newton used it to demonstrate the accuracy of the inverse 
square law in the solar system; notice however that the theorem holds for orbits subject 



to any central force F. Chandrasekhar^ gives an elegant discussion but see Lynden-Bell 
& Lynden-Bell^ for the true shapes of Newton's orbits in uniformly rotating axes. 

Hereafter we specialise to an inverse square law supplemented by an inverse cubic 
force so the potentials considered take the form 

where fi may be thought of as GM and 7^ is a constant which in galactic orbits is 
normally negative though in Hartree Fock atoms it is positive. The equation of motion 
in fixed axes is 

dh/df = -ifxr^'^ + Kr-^)r 

where hats denote unit vectors. 

Comparison with Newton's theorem on revolving orbits suggests that we think of 
the angular momentum h as n~^h and K as (n"^ — l)h'^ = {1 — n'^)h'^ . If we then go to 
axes rotating non-uniformly with angular velocity {n^^ — l)h/r^, the orbital equations 
relative to the rotating axes must reduce to those under the action of the inverse square 
law alone. In practice it is h and K that are known, so in terms of those n^ = 1 — Kh~^. 
Should K be greater than h"^ we would be in trouble, but under such circumstances 
the inverse cube attraction overcomes the centrifugal repulsion so orbits spiral into the 



origin. Otherwise in the non-uniformly rotating axes with f2 = (1 — Vl — Kh 2)h/r^ 
we recover writing a dot for time derivatives relative to the rotating axes 

i^ = — /ir^^f . (1) 

We have indeed laboriously checked that with this f2 the 2 f2 x r+ il x r + f2 x (17 x r) 
terms cancel out with the —Kr/r'^ force term. 

The orbit seen in the rotating axes 

From (dJ 

r X r = h 



a constant, and 



r X h = —/if X ( f X - j = yuf , 



thus 

f X h = /i(f + e) , (2) 

where e is a constant vector fixed in these rotating axes. Since both the other terms are 
perpendicular to h, e must be too, so it hes in the plane of the orbit. Dotting with f//i 
we find setting / = /i^//i, 

l/r = (1 + e.f) = 1 + ecos0 . (3) 

This is the equation of a conic of eccentricity e and semi-latus-rectum / and is the 
angle in the rotating axes measured from perihelion. Thus e points toward perihelion 
in the rotating axes and it will inevitably rotate relative to fixed axes. Since the axes 
rotate at the rate {n~^ — l)h/r'^ = (1 — n)(f) we find 

e = e cos[(l — n)(f)] + h x e sin[(l — n)0] (4) 

where e is in the fixed absolute direction to the perihelion at which = and its 
magnitude is the eccentricity |e|. For fixed axes this comes from Hamilton^. For the 
contributions of Laplace, Runge & Lenz see Goldstein^. 

The Orbit seen from fixed axes 

The equation of the orbit Q is put into fixed axes by writing = n0, / = n'^h'^/fj, 

l/r = 1 + 6 cos n0 . (5) 

The velocity in fixed axes will be 

V = f + f2 X r 

where from Q 

r = ^h X (r + e) . 
nh 



Thus using (Q for e and il = (1 — n)(f)h = (1 — n)h/r^ we find our expression for v 



V = — - < h X r 

nh 



l + -{n-'-l) 



+ h X e cos[(l — n)(j)] — e sin[(l — n)(f)] > (6) 



or using (jSJ to eliminate l/r in favour of 



nh ^- 



nh 



X r 



n ^ + (n ^ — l)ecos(m 



h X ecos[(l — n)(t)\ — esin[(l — n) 



(7) 



To get the radial velocity along any line of sight 1 one merely uses l.v and then 
corrects for the Sun's motion. 



The eccentricity vector constant of the motion 

Equations ® or ((7j) give us the velocity in terms of a conserved constant vector e 
pointing toward the initial perihelion. To get e itself we need to invert this equation so 
e is expressed as a function of v etc. To do this we note that the final two terms in the 
{} in (jH)) yield e if we first cross multiply them by xhcos[(l — n)0] and add the result 
to the same two terms x(— sin[(l — n)(f)\). But from ® those final two terms are equal 

to 

nh - r I . - 

1 



nh 

— V + r X h 



l + -{n 
r 



-1 



hence both e^ is the square of this, which yields e^ = 1 + 2//i ^e where e 
and 



K 
2r^ 
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cos\\\—n)(p\ — < — v+rxn 
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l + -{n-^-l) 
r 



sin[(l— n) 



which gives the new vector constant of the motion. We remind the reader that = 
cos~^(e.f) and n = \/l — Kh~'^. We may eliminate l/r in terms of and e by using 
(0). Notice that (jH} is actually an implicit equation because is not known until e is 
known. However by taking (x, y) coordinates in the plane of the motion and measuring 
a new 0' from the x axis we may write = 0' — 0o where 0o is the azimuth of e. With 
0' V h r 72 / all known it is then possible to find a 0o from the components of (jEJ/e. 



So the implicit equation may be solved for e. There will be very many solutions unless 
we restrict (po to be in the range — vr/ra to +71/71 and (f)' in the range — vr to +7r. I gave 
a discussion earlier^ of the type of forces that leave the magnitude of the eccentricity 
unchanged but slew its direction. 

Approximating orbits in other potentials 

Suppose we are given a potential ipijr) and an orbit within it is defined by its 
pericentre at rp and its apocentre at r^. Then since r vanishes at these points the 
angular momentum of the orbit is given by 

^i){rp) -ipira) 



ly Z. nr» Zi 

p a 



h' = 

and the energy of the orbit is given by 

The angle between perihelion and aphelion is then writing r = u~ 

-1/2 






du 



This may be compared with the angle given by the orbit (0) which is vr/n. 

Hence we may define the n of the approximating orbit by ra = 7r/<l>. 

We shall make the angular momenta of the two orbits equal and define the eccen- 
tricity by 

e = 

Ta + Tp 

which is also in conformity with ^. 

Evidently the K of our approximating potential is already known via n since K = 
h'^{l — n^). To specify the approximating potential we still need fi which we fix by 
making the perihelion distances equal, which, via the eccentricity, implies the aphelion 
distances are equal and hence 

^,= {h'-K){r;' + r-')/2. 
5 



Thus we have an approximation scheme giving n, K, h, e, /i for any orbit specified by r^ 
and Tq in any known potential ^(r). We expect the e defined by (jHJ to be approximately 
constant for such orbits. 

The Eccentricity Vector for the Harmonic Oscillator 

Newton^ explained the connection between the Keplerian ellipse and that generated 
by the two dimensional harmonic oscillator. Here we follow Chandrasekhar's preferred 
path^. Set X + iy = z in the plane of the orbit. The equation of the harmonic oscillator 
is then ^ + u'^z = and its energy is E = ^{^^ + uP''zz). Its angular momentum is 
^ ~ h.^%. ~ ^'^)' N*^^ consider the mapping of the complex plane Z = z^ . Following 
Newton we ask whether the mapped path considered with a new time r(t) can be an 
orbit under a new central force. Evidently the angular momentum of the Z orbit is 

1 (—dZ dZ \ 1 ,2[-dz dz \ dt ^ 

— z— —Z = -\z\ z- —z = 2-—\z\h 

2i \ dr dr J i \ dr dr J dr 

so if the angular momenta are to be equal then ^ = g^F^t' Chandrasekar (in error!) 
omits the factor 2. Now using the z equation of motion 

d'^Z 1 d (2dz\ 1 (d-zdz o_\ s eZ 

— ijj zz 



dr^ A\z'^\ dt \-z dt J 2zz^ \dt dt J zz^ \Z\^ 

but the final expression shows us Z is a motion under an inverse square law with a 
force constant fi = GM = e where e is the energy of the simple harmonic orbit. Thus 
under the mapping the simple harmonic centred ellipse becomes the Kepler eccentric 
ellipse. However the latter has a conserved Hamilton eccentricity vector so what does 
that vector become under the inverse transformation from Kepler's to Hooke's ellipse? 
In the notation of this section vectors in the plane of the motion are complex numbers. 
Since d^Z/dr"^ = —fiZ/\Z\^ we follow our well trodden path and multiply by h and 
integrate to find 

hdZ/dT = iii{Z/\Z\+e) 

where e is the (complex) constant of integration. To transform this into the z plane we 

6 



write d/dr = gUp^, Z = z^ and jjl = e. So 

-[---n + -]=e. (9) 

\ e z dt z j 

I had to differentiate this extraordinary expression and show the resuh to be zero 
before I beheved that it was indeed a constant of the motion! In the original Z space e 
pointed toward perihehon. After the transformation Z = z^, e is still the same complex 
number so is unchanged but now the perihelion in z space will be at half the angle to 
the real axis. This suggests that we should be considering a new vector e^ pointing to 
perihelion and with the property that e = ei. This e^ will then be a constant of the 
motion too and furthermore the ± ambiguity in its definition reflects the fact that the 
Hooke ellipse has two perihelia in opposite directions. There is an intrinsic difficulty in 
the transformation that we have to place a cut in the complex Z plane to define \fZ 
properly. It we place this cut arbitrarily then the direction of that cut intrudes into the 
resultant formulae. It is much more sensible to take the cut to be defined physically. 
Taking the cut along the real axis and that toward perihelion along e has advantages. 
Then both e and Cp are real. Rewriting our equation to give us the velocity we have 

dz ie _. 

we may now rewrite this in vector form 

— = - {h X [r + ep.rep + Gp x (gp x r)]} 

which gives the velocity in terms of the vector Gp that points to perihelion and the 
position vector r in the orbital plane orthogonal to h. To find the magnitude of e 
and Cp we return to equation Q. The general solution to the harmonic equation is 
z = pe"^*+ge~*'^* where p and q are complex numbers. In terms oip and q, h = uj{pp~qq) 
and e = a;^(pp+ qq). Putting these expressions into equation (jH)) we find 

—2pq 

e = -^ 

pp + qq 



choosing the real axis along e means that p = Pe^^ and q = —Qe ^^ with P and Q real 
and positive, then z = {P — Q) cos(x + ^t) + i{P + Q) sin(x + ojt) 



2PQ a' - ¥ e 



P2 + g2 a2 + ^2 2 - e2 ' 
where a, h are the semi axes and e* is the 'eccentricity' of the centred Hooke ellipse. 
Thus the magnitude of Cp is given by 



Cp 



2 -el 



Of course e^ could be generalised to cases in which the linear Hooke law is supple- 
mented by an inverse cube repulsion, by following the method given in this paper. 
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